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We derive exact relations between the Potts model partition function, or equivalently the Tutte 
polynomial, for a network (graph) G and a network obtained from G by (i) by replacing each 
edge (i.e., bond) of G by two or more edges joining the same vertices, and (ii) by inserting one 
or more degree-2 vertices on edges of G. These processes are called edge and vertex inflation, 
respectively. The physical effects of these edge and vertex inflations are discussed. We also present 
exact calculations of these polynomials for families of networks obtained via the operation (ii) on 
a subset of the bonds of the network. Applications of these results include calculations of some 
network reliability polynomials. In addition, we evaluate our results to calculate various quantities 
of structural interest such as numbers of spanning trees, etc., and to determine their asymptotic 
behavior for large networks. 
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I. INTRODUCTION 

The Potts model has long been valuable as a system 
exhibiting many-body cooperative phenomena On 
a lattice, or, more generally, on a network (i.e., graph) 
G, at temperature T, the partition function for this 
model (in zero external field) is Z = J2{a } e^*^^, where 
/3 = l/^ksT), the Hamiltonian H = —JJ2e '^ffi^ji J is 
the spin-spin exchange constant, i and j denote vertices 
(sites) on G, is the edge (bond) connecting them, and 
tJi are classical spins taking values in the set 
We use the notation K = /3J and v = — 1. Thus, for 
the Potts ferromagnet (FM, J > 0) and antiferromagnet 
(AFM, J < 0), the physical ranges of v are v > and 
— l<v<0, respectively. In general, a graph G = {V, E) 
is defined by its vertex set, V, and its edge set, E. We de- 
note the number of vertices of G as n = n(G) = \V\ and 
the number of edges of G as e(G) = Z is a. polyno- 
mial in q and v, as will be evident from Eq. (|2.ip below. 
The Potts model partition function Z is equivalent to a 
function of considerable interest in modern mathematical 
graph theory, namely the Tutte polynomial '^j-fJ], and 
we shall therefore often refer to this object in a unified 
manner as the Potts/Tutte polynomial. The partition 
function of the zero-temperature Potts antiferromagnet 
is identical to another function of longstanding interest 
in graph theory, namely the chromatic polynomial, which 
counts the number of ways of assigning q colors to the 
vertices of G subject to the condition that no two adja- 
cent vertices have the same color. These are called proper 
g-colorings of G. An important property of the Potts an- 
tiferromagnet is that for sufficiently large q, it exhibits 
nonzero entropy per site at zero temperature and is thus 
an exception to the third law of thermodynamics @ . 

An interesting problem involving both statistical me- 
chanics and mathematical graph theory is to relate the 
Potts/Tutte polynomial calculated for a graph G with 



the corresponding polynomial calculated for a graph G 
which is obtained by a specified modification of G. In 
this paper we will first present a general solution to this 
problem for two classes of G's, namely those obtained 
(i) by replacing each edge of G by two or more edges 
joining the same vertices, and (ii) by inserting one or 
more degree-2 vertices on each edge of G (where the de- 
gree, tvu - , of a vertex Vi is defined as the number of edges 
that connect to it). We denote these operations as edge 
and vertex inflations of G, respectively. In the litera- 
ture on mathematical graph theory, a vertex inflation 
of G is also called a homeomorphic inflation of G, and 
the reverse procedure, of removing one or more degree-2 
vertices from edges of G, is called a homeomorphic re- 
duction of G. In the remaining sections of our paper 
we will present results for Potts/Tutte polynomials of 
certain families of graphs obtained by vertex inflations 
for a subset of edges. These families include longitudi- 
nal vertex inflations of free, cyclic, and Mobius ladder 
graphs of arbitrarily great length and a family that we 
call hammock graphs. These results extend our previ- 
ous work in two ways: (a) as generalizations to the full 
Potts/Tutte polynomials of our earlier calculations with 
S.-H. Tsai of the chromatic polynomials for homeomor- 
phic inflations of these families of graphs in Refs. 0-0 
and (b) as homeomorphic inflations of our previous cal- 
culation of the Potts/Tutte polynomials for ladder strips 
inRef. 

There are several motivations for this work. For an ar- 
bitrary graph G, the calculation of the Potts/Tutte poly- 
nomial involves a number of computational steps and a 
corresponding time that grow exponentially rapidly with 
n(G) and e(G) (e.g., 11, 12]). Furthermore, there is no 
known exact closed-form solution for Z for arbitrary q 
and temperature T on the thermodynamic limit of lat- 
tice graphs of dimensionality d > 2. Hence, it is of funda- 
mental value to carry out exact analytic calculations of 
Potts/Tutte polynomials on various families of graphs. 
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such as lattice strip graphs and modifications thereof. 
Furthermore, special cases of the Potts/Tutte polynomial 
are of considerable interest in their own right. We have 
noted the importance of the zero-temperature Potts anti- 
ferromagnct (chromatic polynomial). Thus, another mo- 
tivation for our work, as embodied in point (a) above, is 
to generalize our previous calculations with S.-H. Tsai [7|- 
to the broader context of the finite-temperature Potts 
antiferromagnet and the Potts ferromagnet (the latter 
also at arbitrary temperature). A central motivation is 
to investigate the effects of edge and vertex inflations of 
a graph on the Potts/Tutte polynomial of that graph. 
This involves the generalization in point (b) above. One 
particular value of our present work is its demonstration 
of the usefulness of relating calculations in the context 
of the Potts model to corresponding calculations for the 
Tutte polynomial and vice versa. Finally, special cases 
of our results yield various quantities of interest such as 
reliability polynomials, numbers of spanning trees, etc. 
for these families of graphs. 



II. SOME GENERAL BACKGROUND 

In this section we briefiy discuss some necessary back- 
ground that will be used for our calculations. There is 
a useful relation that expresses the Potts model parti- 
tion function Z{G,q,v) as a sum of contributions from 
spanning subgraphs of G. Here, a spanning subgraph 
G' ~ {V,E') has the same vertex set as G and a subset 
of the edge set E, E' C E. This relation is [11] 



Z{G,q,v) 



^ ^fcc(G')^e(G') 
G'CG 



(2.1) 



where kf,{G') denotes the number of connected compo- 
nents in G' . This formula shows that Z{G^ q, v) is a poly- 
nomial in q and v. For the Potts ferromagnet, Eq. (12. ip 
enables one to generalize q from the non- negative integers 
to the non-negative real numbers while keeping Z(G, q, v) 
positive and hence maintaining a Gibbs measure. 

The Tutte polynomial of a graph G, denoted 
T{G,x,y), is defined by Q-jl 

T(G,x,2/)= ^ (.T~l)'==(°')-'^'=(^)(y-l)^(^') , (2.2) 

G'CG 

where c(G") is the number of linearly independent cycles 
(circuits) in G". Let us define 



q 

X = 1 ^ — , y = e =w + l 



(2.3) 



so that (7 = (x— l)(j/— 1). The relation between Z(G, q, u) 
and r(G, x, y) follows directly from Eqs. ((2?T|) and (|2T2)) 
and is 

Z(G, g, v)^{x- 1)'=^(G) [y - 1)"(G)T(G, x, y) . (2.4) 

With no loss of generality, we will restrict here to con- 
nected graphs G, so fcc(G) = 1. 



For a graph G, let us denote G— e as the graph obtained 
by deleting the edge e and G/e as the graph obtained by 
deleting the edge e and identifying the two vertices that 
were connected by this edge of G. This operation is called 
a contraction of G on e. From Eq. (|2.1I) . it follows that 
Z{G, q, v) satisfies the deletion-contraction relation 

Z{G, q, v) ^ Z{G - e, q, v) + vZ{G/e, q, v) . (2.5) 

An analogous deletion-contraction relation holds for 
r(G,x,y). 



III. TYPES OF INFLATIONS OF A GRAPH 

Here we discuss in greater detail the edge and vertex 
inflation of a graph G. As part of our study, we will con- 
sider the infinite-length limit of lattice strip graphs G of 
finite width, and also the thermodynamic limit of lattice 
graphs G of dimensionality d > 2. In these cases, we 
will often use the notation {G} to indicate these limits. 
With no loss of generality, we begin by assuming that G 
has no multiple edges. We define an edge infiation of G 
to be a graph obtained by replacing one or more of the 
edges of G by multiple edges joining the same vertices. 
Clearly, this leaves the number of vertices n(G) invariant. 
In particular, it is natural to define a uniform ^-fold edge 
inflation of G as the graph ff(G) obtained by replacing 
each edge of G by £ edges joining the same two vertices. 
Thus, the number of edges of £e{G) is e{£i{G)) = ie{G). 
A K-regular graph is a graph with the property that all 
vertices have the same deg y eV. Al- 

though a general graph is not K-regular, one can define 
an average or effective vertex degree Hgf / as 



y 



n{G) 



(3.1) 



Clearly, a uniform ^-fold inflation of all of the edges of G 
multiplies K^ff by i. A remark is in order here concern- 
ing loops. A loop is defined as an edge joining a vertex 
back to itself. This would not occur in the statistical 
mechanical framework, since it would mean a spin ct; in- 
teracting with itself, and hence we will usually assume 
that G is a loopless graph. 

We define a vertex (i.e., homeomorphic) inflation of G 
to be a graph obtained by inserting one or more degree-2 
vertices on one or several edges of G [l3|. This leaves 
the number of (linearly independent) circuits in G, c(G), 
invariant. Homeomorphic inflations and reductions of 
graphs have been of interest in the study of chromatic 
and Tutte polynomials and have also been 

studied in the context of "decorated" spin models [23]. 
As with edge inflations, it is natural to define a uniform 
^-fold vertex infiation Vi{G) as the graph obtained by 
inserting £ degree-2 vertices on each edge of G. For the 
special case where G is a section of a regular lattice, it is 
also natural to consider edge or homeomorphic inflation 
of edges forming a subset of lattice vectors. For example. 
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for a d-dimensional Euclidean lattice E*^, one can consider 
the case in which the edge or vertex inflation is performed 
on the edges along one or more lattice vectors Cj, where 
j takes values in a subset of {1, d}. 



IV. RELATIONS FOR UNIFORM EDGE AND 
VERTEX INFLATIONS 

A. Edge Inflation 

The effect of a uniform edge inflation on the Potts 
model partition function can be determined from an anal- 
ysis of the Potts Hamiltonian. As we edge-expand G to 
£e{G), i.e., replace each edge by £ edges joining the same 
vertices, we induce the change in the Hamiltonian 



n 



i.e., J — >■ £J, and hence 

y Ve.e = y'^ 



(4.1) 



(4.2) 



where the subscripts e, £ refer to the £-fold edge inflation. 
Equivalently, 



V Ve,i = {v + lY -1 , 



(4.3) 



where v^^ g = y^^^ — 1. This proves the following relation 
connecting Z on G with Z on £e{G): 



Z{£e{G),q,v)^Z{G,q,Ve^i) 



(4.4) 



We next determine the effect of this edge inflation on 
the Tutte polynomial. Given the transformation (|4.2p 
and the fact that q does not change, so that 

q^ix- l){y - 1) = (xej - l){ye,i - 1) , (4.5) 

we have 

' x-1 



^ , ix~l){y-l) 

Xe,e = IH -f — = 1 

2/ - 1 



(4.6) 



Combining these transformations of variables with Eq. 
(12. 4p . we derive the relation connecting the Tutte poly- 
nomials of G and of £e{G), namely 



T{£^{G),x, y) = T{G, Xe,e, ye,e) 



(4.7) 



where x^^i and yej were given in Eqs. (14. 6p and (j4.2l) . 
Note that the effect of the edge inflation on the Potts 
partition function is simpler than the effect on the Tutte 
polynomial, since in the former case, only one of its vari- 
ables is modified, namely, v — > Ve,e, whereas in the latter 
case, both of its variables are modified, as a; — )■ Xe,i and 
y ye/. 

Although our main focus is on the Potts model with 
spin-spin exchange constants J that are independent of 



the edges Cij, we note parenthetically that one can con- 
sider a generalization in which the Jij are different for 
each edge of G, e^. From the original 1944 Onsager 
solution of the two-dimensional Ising model {^S*] , a large 
number of studies of spin models have dealt with the gen- 
eral case with different spin-spin exchange constants for 
different lattice directions. Studies of spin models with 
spin-spin exchange constants Jij that can be different (in 
magnitude and sign) for each edge, were motivated by 
early work on spin glasses [13] . In this case, where Jij 's 
depend on , the transformation of the Hamiltonian be- 
comes 



(4.8) 



Hence, defining Kij = fiJij and Vij = e^*^ — 1, we have 

Vij Vij,e,e = [vij -f 1)^ - 1 . (4.9) 

Denoting {u} as the set of w^'s, we then have 

Z(f,(G), g, {v}) = Z(G, q, {v,A) . (4.10) 

Since J < for the Potts antiferromagnet, it follows 
that, as T — >■ 0, — >■ — oo (and v — >■ —1). Hence, in this 
limit, the only contributions to the partition function are 
from spin configurations in which adjacent spins have dif- 
ferent values. The resultant T = Potts AFM partition 
function is therefore precisely the chromatic polynomial 
P(G, q) of the graph G counting the number of proper 
g-colorings of G: 



Z(G,g,-l)=P(G,g) . 
This is equivalent, via Eq. (|2.4I) . to the relation 



(4.11) 



P(G, q) = (-g)'==(«)(-l)"(G)T(G, 1 - g, 0) . (4.12) 

The minimum number of colors necessary for a proper 
g-coloring of G is the chromatic number x(G). For 
q > x(G), P{G,q) grows exponentially with n, leading 
to a ground state degeneracy per vertex, W{{G},q) > 1, 
where T^({G},g) = lim„^oo P(G, g)i/". The ground 
state entropy per vertex of the Potts model on |G} is 
So{{G},q) = kBHW{{G},q)]. In Refs. Q, [3, 
we applied our exact results on chromatic polynomials to 
study the phenomenon of nonzero ground state entropy 
per site in Potts antiferromagnets. These exact results 
complement other approaches to studying 5*0, such as 
rigorous bounds, series, and Monte Carlo measurements 

It is clear from the definition of the chromatic poly- 
nomial that P(G, q) does not change if one replaces any 
edge of G by two or more edges joining the same vertices. 
In particular, for the case of an ^-fold uniform edge in- 
flation of G, 



Pi£iiG),q)^PiG,q) 



(4.13) 



This is also clear analytically from Eqs. (|4.3p and (14.41) . 
since the condition that v = —1 implies that v^.e = — 1. 
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B. Vertex Inflation 

To analyze the effect of a vertex inflation of a graph G, 
it is again convenient to start with the Potts model for- 
mulation. For simplicity, we assume here that G does not 
have any multiple edges; it is straightforward to extend 
our calculation to the case of multiple edges. We use the 
fact that in the basic expression Z = J2{a } ^ '^^^ 
can perform the summations over the spins that are lo- 
cated at these degree-2 vertices. Let us consider an edge, 
Cij, and insert a degree-2 vertex Va (and its associated 
spin, a a) on this edge. Then 



(Ti.(7a-(7 



(4.14) 

Now, carrying out the summation over aa, we find that 
if CTi = (Tj , then the result is q + 2v + v'^, while if ai ^ aj , 
then the result is 5 + 2v. Therefore, 

= {q + 2v)Y,{^ + Vvs5„^.,) , (4.15) 



where 



Vv.l 



q + 2v 



(4.16) 



and the subscript v, 1 refers to the insertion of one addi- 
tional vertex on the edges. Performing this summation 
for each edge, we derive the relation 

Z(Vi(G),<z,t;) = (9 + 2«)^(«)Z(G,g,«,„.i) . (4.17) 

This operation can be performed iteratively i times, 
thereby giving a relation between Z on Vi{G) and Z on 
G. For example, for £ = 2, we have 



Vv,2 



q + 2vy^i 



{q + 2v){q^ + 2qv + 2v'^ 



(4.18) 



and so forth for higher values of We can thus relate Z 
on Vi(G) to Z on G. 

A basic problem in graph theory is the enumeration 
of discretized flows on the edges of a (connected) G that 
satisfy flow conservation at each vertex, i.e. for which 
there are no sources or sinks. One arbitrarily chooses 
a direction for each edge of G and assigns a discretized 
flow value to it. The value zero is excluded, since it is 
equivalent to the edge being absent from G; henceforth, 
we take a q-flow to mean implicitly a nowhere-zero q-flow. 



The flow on each edge can thus take on any oiq—1 values 
modulo q. The flow or current conservation condition is 
that the flows into any vertex must be equal, mod g, to 
the flows outward from this vertex. These are called q- 
flows on G, and the number of these is given by the flow 
polynomial, F{G,q). This is a special case of the Tutte 
polynomial for a; = and y = 1 ~ q or equivalently, in 
terms of Potts model variables, v — —q: 



F{G,q) = {-l)<^^TiG,0,l-q) 



(4.19) 



As is clear from the definition of the flow polynomial, 
adding or removing a degree-2 vertex from an edge of G 
does not change the number of allowed q-flows on G, so 



F{V,{G),q)^F{G,q) . 



(4.20) 



As is evident from Eq. (I4.16|) . the condition v — —q 
implies that Vy^i = —q and hence, more generally, that 
Vv,e — —q for arbitrary £ > 1. Recall that the number of 
cycles in G is unchanged by this homeomorphic inflation: 
c{Vi{G))^c{G). 



V. PHYSICAL EFFECTS OF EDGE AND 
VERTEX INFLATION 

A. General 

An important question concerns the effect of edge and 
vertex inflations on the physical properties of the Potts 
model. First, one may investigate these effects for the 
q-state Potts model on (the thermodynamic limit of) a 
regular lattice graph of dimension d > 2, where the fer- 
romagnetic version of the model has a flnite-temperature 
order-disorder phase transition, and, depending on the 
lattice type and the value of q, the Potts antiferromag- 
net may have a finite-temperature order-disorder phase 
transition. In particular, one can study the effect on the 
phase transition (pt) temperature Tpt for lattices where 
the equation for this quantity is known exactly. How- 
ever, aside from the q = 2 Ising case on two-dimensional 
lattices, there is no known exact closed-form solution 
for the free energy of the q-state Potts model at arbi- 
trary temperature on these lattices with d > 2. There 
is thus also some interest in studying the effect of edge 
and vertex inflation for inflnite-length limits of quasi-one- 
dimensional lattice strips, where one can obtain exact 
closed-form expressions for the free energy for arbitrary 
q and T. Of course, a spin model with short-ranged in- 
teractions, such as the Potts model, does not have any 
flnite-temperature phase transition on an inflnite-length 
quasi-one-dimensional lattice strip. Nevertheless, one in- 
teresting application of exact solutions for the free energy 
and thermodynamic quantities on these strips is that one 
can study their dependence (and the dependence of the 
associated Tutte polynomials) on graphical properties, in 
particular, on the edge or vertex inflation. 

In accordance with our notation above for the dimen- 
sionless inverse temperature K = J/{kBT), we define 
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pt 



J /{ksTpt)- We also introduce some notation transition in the presence of greater thermal fluctua- 



that we will use below. We define the shifted values 
of Tpt due to a uniform ^-fold edge inflation (symbol- 
ized by the subscript e,t} and a uniform £-fold vertex 
inflation (symbolized by the subscript by Tpt^e,i 

and Tpt.v,i- We thus denote Kpt.e,e = J/iksTpt^ej) and 
Kpt^vi = J/{kBTpt.v,e)- 



B. Effect on Tpt Due to Edge Inffation 

A general result is that for the Potts model on (an infi- 
nite) regular lattice graph {G} with dimensionality d> 2 
that has a finite-temperature phase transition (of cither 
ferromagnetic or antiferromagnetic type, and either first 
or second order) at a temperature Tpt, a uniform ^-fold 
inflation of all edges of {G} has the effect of multiplying 
Tpt by the factor £, i.e., 



T 



iT 



pt 



K, 



pt.e,i 



pt 



(5.1) 



Analytically, this follows because the uniform i-io\d edge 
inflation of {G} changes J to £J and Tpt is proportional to 
J. Physically, it follows because replacing J by £J with 
£ > 2 strengthens the spin-spin interaction and hence 
makes possible the onset of long-range magnetic order 
in the presence of greater thermal fluctuations, i.e., at a 
higher temperature. 

As an example, consider the Potts ferromagnet on the 
(thermodynamic limit of the) square lattice, {sq} = E^. 
Denoting the spin-spin exchange constants in the two lat- 
tice directions e^, i — 1,2, as Ji, with Ki = /3Ji and 
Vi = e^' — 1, we recall the well-known equation for the 
phase transition temperature, namely [l|. 



V1V2 = q ■ 



(5.2) 



Let us initially assume Ki — K2 = K and thus vi — 
V2 = w, so Eq. (j5.2l) becomes = q. The dimensionless 
inverse phase transition temperature Kpt is then given 
by 



Kpt = ln(l + 



(5.3) 



Now let us carry out an £-fold edge inflation on the lattice. 
Denoting the resultant inverse phase transition temper- 
ature in an obvious notation as Kpt.e,e, we have 



pt,e,t = J ln(l + Vq) = ^ . 



(5.4) 



One can also consider the effect of an edge inflation 
on all edges along a subset of lattice directions. The 
effect is simplest for the ferromagnetic case, since this 
does not involves competing interactions or frustration. 
This edge inflation along a subset of lattice directions 
strengthens the net spin-spin interaction and therefore 
makes possible the ordering associated with the phase 



tions. Let us denote Tpt^es,i as the shifted phase tran- 
sition temperature after an £-fold edge inflation along 
a subset s of the lattice directions, and similarly de- 
note Kpt^es,i = J / {kBTpt,es,i)- Thcn the reasoning above 
yields the inequality Tpt^es,e > Tpt for £> 2. As an exam- 
ple, we again consider the (infinite) square lattice {sq} 
and perform an edge inflation with £ = 2 for edges in 
one of the two lattice directions, say 62, so that Ki = K, 
K2 = 2K. The equation for the shifted inverse phase 
transition temperature is given by v'^{v + 2) — q, with 
physical solution 



Kpt^e2,2 — In 



3l 



(5.5) 



where the subscript e2 means inflation along edges along 
62 and 



A 



1 r 



27q - 16 + 3V3g(27q^^32) 



(5.1 



The resultant Kpt^e2,2 < Kpt, i.e., Tpt^e2a > Tpt, in agree- 
ment with the general argument given above. For exam- 
ple, for q^2, Kpt = \n{l + V2 ) ~ 0.88137, while i^pt,e2,2 
is given by 



SO that K 



(19-F3V33)i/^-F4(19-H3V33 )"i/3-M , 

(5.7) 



pt,e2,2 



0.60938. 



C. Effect on Tpt Due to Vertex Inffation 

One can also deduce a general result for the effect of a 
uniform £-fold vertex (i.e., homemorphic) inflation of the 
thermodynamic limit of a lattice graph {G} with dimen- 
sionality d > 2. A generic feature of the phase transition 
temperature Tpt of a ferromagnetic spin model (above its 
lower critical dimensionality, so that this temperature is 
flnitc) is that, other things being equal, Tpt increases as 
a function of the vertex degree, i.e., coordination num- 
ber, of the lattice. This feature is observed in approx- 
imate determinations of Tpt from high-temperature and 
low-temperature series expansions, Monte Carlo simula- 
tions, mcan-fleld approximations, and, where available, 
exact solutions for Tpt- It is understood physically as a 
consequence of the fact that increasing the coordination 
number increases the effect of the spin-spin interactions, 
so that the ordering associated with the phase transition 
can occur in the presence of greater thermal fluctuations. 
On an (inflnite) line, with coordination number k = 2, a 
spin model with short-range interactions does not have 
a finite-temperature phase transition, so the lattices of 
interest in this section are lattices with d > 2, which 
necessarily have coordination number k > 3 (where this 
minimum value, k = 3, is realized for the honeycomb 
lattice). If one starts with a lattice graph {G} with coor- 
dination number k > 3, then a uniform vertex infiation. 
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which consists of the addition of i degree-2 vertices on 
each edge of {G}, reduces the effective vertex degree, 

Let us, for technical simpHcity, consider the thermo- 
dynamic hmit {G} to be reached as the n{G) -> oo 
hmit of a regular lattice graph G with periodic bound- 
ary conditions, and uniform vertex degree k. Then, with 
e{G) = {n/2)n{G), the number of vertices and edges of 
Vi{G) are 



and 



so that 



n(V,(G))=(l + ^)n(G) 



e(V,(G)) = (l+^)e(G) 



«e//(V,(G)) = 



1 + ^ 



(5.8) 



(5.9) 



(5.10) 



The fact that the vertex inflation reduces Heff if k > 2 
is clear analytically from this result, since 



Ke//(Vf(G)) < K if £ > 1 and k > 2 



(5.11) 



For example, for the square lattice, with k = 4, one has 
K,f}{Vi{sq)) - 8/3 - 2.666.., n,f}{V2{sq)) = 12/5 = 
2.4, Ke//(V3(s(7)) = 16/7 ~ 2.286, and so forth, with an 
approach to 2 from above as £ — ^ cxd. Indeed, in general, 
for an arbitrary regular lattice graph G with coordination 
number k. 



lim K 



For i >> 1, one has the Taylor series expansion 



= 2 



2 



(5.12) 



(5.13) 



Hence, for a Potts ferromagnet on a lattice graph {G} 
with d > 2, a uniform ^-fold vertex inflation of {G} with 
i > 1 leads to a decrease in the phase transition tem- 
perature. The same conclusion holds if one performs a 
vertex inflation on all edges along a subset of the lattice 
directions. 

We illustrate the effect of vertex inflation for the q- 
state Potts ferromagnet on the square lattice. Let us 
perform a uniform vertex inflation with ^ = 1 on all edges 
of the lattice, i.e., add one degree-2 vertex to each edge 
of this lattice. Then, combining Eqs. (|4.3|) and (|5.2p . we 
find that the equation for the phase transition temper- 
ature is Vv^i = y/q, where Vy^i was given in Eq. (|4.16|) . 
The solution for the inverse phase transition temperature 



= ln 



(5.14) 



Comparing this with the inverse critical temperature for 
the original {G}, Kpt = ln(l-|-^ ), we see that Kpt^v,i > 
Kpt, in agreement with our general argument above. 

The situation is more complicated with the Potts an- 
tiferromagnet; we will show that vertex inflation can ei- 
ther lower or raise a phase transition (critical) temper- 
ature, depending on the value of q and the lattice type. 
First, consider the q = 2 (Ising) Potts antiferromagnet 
on a bipartite lattice {G}. The bipartite property of 
{G} means that it can be expressed as the union of an 
even and an odd sublattice, {G} — {Gi} U {G2}, with 
the property that each vertex in Gi has, as its only ad- 
jacent vertices, members of the vertex set of {G2} and 
vice versa. There is a well-known isomorphism that maps 
the Ising antiferromagnet on {G} to an Ising ferromagnet 
on {G}, namely the simultaneous replacement J — >■ — J 
and (T^j^ — > — (Jtjj , with a^^ unchanged, where here vi and 
V2 denote vertices in Gi and G2, respectively. Thus, in 
this case, the effect of an £-fold vertex inflation on all 
edges is the same for the antiferromagnetic case as for 
the ferromagnetic case discussed above, namely that it 
reduces Tpt. However, we next prove that vertex infla- 
tion can also have the opposite effect, of increasing a 
critical temperature. For this purpose, let us consider 
the q = 2 Potts (Ising) antiferromagnet on the infinite 
triangular lattice {tri} (with equal negative spin-spin ex- 
change constants in each of the three lattice directions, 
Ji = J2 = J3 = J < 0). This antiferromagnet is frus- 
trated, and is only critical at T = \2^. Now let us per- 
form a uniform £-fold vertex inflation on all edges, with 
i — 2k + 1 odd, thereby obtaining {V2k+i{tri)}. This 
lattice, {V2fe+i(iri)}, is bipartite, in contrast with the 
triangular lattice itself. Because of this, the Ising antifer- 
romagnet is not frustrated on {V2fe+i(iH)}, and therefore 
one can apply a standard Peierls-type argument to infer 
that it has a finite-temperature symmetry-breaking phase 
transition. Indeed, because {V2k+i{tri)} is bipartite and 
because of the isomorphism mentioned above, the Ising 
ferromagnet and antiferromagnet can be mapped to each 
other, and have their respective phase transitions at the 
same Tpt. Thus, as this example shows, in contrast with 
the situation for the Potts ferromagnet, vertex inflation 
for the Potts antiferromagnet may actually raise a criti- 
cal or phase transition temperature rather than lowering 
it, depending on q and the lattice type. 



D. Invariance of Universality Class Under Edge 
Inflations 



We recall that on two-dimensional lattices the (zero- 
field) g-state Potts ferromagnet with g < 4 has a second- 
order phase transition with an associated g-dependent 
universality class and corresponding thermal and mag- 
netic critical exponents that are independent of the lat- 
tice type. It is of interest to study whether vertex or edge 
inflation changes the universality class of this phase tran- 
sition. A general result of renormalization-group analyses 
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of second-order phase transitions is that the universal- 
ity class of a second-order phase transition (in a model 
that is free of complications such as competing interac- 
tions, frustration, and/or quenched disorder) depends on 
the lattice dimensionality, d and the symmetry group of 
the Hamiltonian (e.g., |28]). For the Potts model, the 
symmetry group of the Hamiltonian is the symmetric 
(permutation) group on q indices, Sq. Edge inflations of 
the lattice have no effect on the dimensionality or sym- 
metry group of the Hamiltonian; hence, for the range 
q < 4: where the two-dimensional Potts ferromagnet has 
a second-order phase transition, they do not change the 
universality class of this transition. 



E. Effects of Vertex Inflation on Universality Class 

We again consider the interval g < 4 where the two- 
dimensional Potts ferromagnet has a second-order phase 
transition. As with edge inflation, we note that vertex 
inflation has no effect on the lattice dimensionality or 
symmetry group of the Hamiltonian, and hence, by the 
same argument as before, it does not change the univer- 
sality class of the phase transition. 

As before, the situation is more complicated for the 
Potts antiferromagnet, because, in contrast to the ferro- 
magnet, its properties depend sensitively on the type of 
lattice. To show this, it is convenient to continue with 
the same illustrative example that we used above, namely 
the q = 2 Ising antiferromagnet. We will give two exam- 
ples which exhibit opposite behaviors; in the first, the 
vertex inflation lowers the phase transition temperature 
and makes no change in the universality class. In the sec- 
ond, the vertex inflation raises the critical temperature 
and does change the universality class. The first exam- 
ple uses the Ising antiferromagnet on a bipartite lattice 
graph {G} of dimensionality d> 2, where it has a phase 
transition at a temperature Tpt (with antifcrromagnetic 
long-range order for T < Tpt). Now let us perform a 
uniform £-fold vertex inflation on all of the edges of {G}, 
thereby obtaining {Ve{G)}. By an argument similar to 
the one given above, this lowers Tpt . Since this vertex in- 
flation does not change either the lattice dimensionality 
or the symmetry group 5*2 ~ , it leaves the universality 
class unchanged. 

To show that the opposite can also happen, let us con- 
sider the (isotropic) Ising antiferromagnet on the trian- 
gular lattice. As mentioned before, because of the frus- 
tration, this model has no finite-temperature transition, 
but is critical at T = 0; the spin-spin correlation function 
decays asymptotically like (aocrp) cx r^^/^ cos(27rr/3) for 
large r = jf] [27j. Normally, at a second-order phase 
transition of a spin model on a d-dimensional lattice in 
which the (connected) spin-spin correlation function de- 
cays asymptotically like (aocrf) oc r~'-'^~^^^), one assigns 
the critical exponent rj to this transition. Because of the 
oscillatory nature of the asymptotic decay of the T = 
Ising antiferromagnet on the triangular lattice, this case 



is more complex, but the decay of the envelope curve 
is described by 77 = 1/2. Now, just as we did before, 
let us perform a uniform ^-fold vertex inflation on all 
edges, with £ — 2k + 1 odd, thereby obtaining the bi- 
partite lattice {V2k+i{tri)}. As discussed above, because 
{V2A:+i(tri)} is bipartite, the Ising antiferromagnet is not 
frustrated on it, and, indeed, can be mapped to the Ising 
ferromagnet by the mapping given in the previous sub- 
section. Owing to this, the Ising ferromagnet and antifer- 
romagnet on this lattice have the same phase transition 
temperatures, and, furthermore, the Ising antiferromag- 
net is automatically in the same universality class as the 
Ising ferromagnet, with = 1/4 [2^, [2^. Thus, in this 
case, the vertex inflation does change both the value of 
the critical temperature and the universality class of the 
phase transition. 



VI. LONGITUDINAL HOMEOMORPHIC 
INFLATIONS OF FREE LADDER GRAPH 

In the previous sections we have derived general re- 
lations that connect Z{G,q^v) and Z{G, q,v), where G 
is obtained from G by uniform edge or vertex inflations. 
By Eq. (|2.4|) . these enable one to calculate the equiva- 
lent Tutte polynomial, T(G, x, y). We have also discussed 
the case where edge or vertex inflations are performed on 
all edges along a subset of lattice directions of a lattice 
graph. In the rest of this paper we explore the latter 
type of vertex inflation further. We present exact calcu- 
lations of Potts/Tutte polynomials for a class of vertex 
(i.e., homeomorphic) inflations of a subset of the edges 
of ladder graphs. Our results generalize our calculations 
of the chromatic polynomials for the se g raphs with S.-H. 
Tsai in Refs. [| and (see also i,|T6|). 

We consider the free strip of the ladder graph com- 
prised of m squares, which we denote as Sm- Now we 
perform an ^-fold vertex inflation on all of the longitudi- 
nal edges, with I — k — 2 and fc > 3, i.e., we insert k — 2 
degree-2 vertices on each of these edges. The resultant 
strip graph is denoted St.m- The transverse edges are 
not affected by this operation. The original ladder graph 
itself is S'2,m- The numbers of vertices and edges on Sk.m 
are 

n{Su.,m) = '2{k-l)m + 2 (6.1) 

and 

e{Sk,„,) ^ {2k - l)m + 1 . (6.2) 

Using a systematic iterative application of the deletion- 
contraction property, we calculate the Tutte polynomial, 
T{Sk,rm x,y). One way to express this is in terms of a 
generating function. We use a generating function 

00 

r{Sk,x,y;z) ^ ^T{Sk.m+i:X,y)z"' (6.3) 

m=0 
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of the form 



ai = -yx 



2k-l 



(6.8) 



T{Sk,x,y;z) 



ao + aiz 



(6.4) 



The denominator can be written as 

l + biz + b2Z^ = (1 - A,,„,z)(l - A,,, ,z) . (6.5) 

Recall that for the circuit graph with n vertices, C„, 

T{C,„x,y)^ ^ + (6.6) 
X — 1 ^-^ 

i=i 

where c^^^ = q— l = xy — x^y. We calculate 

ao = T{C2k,x, y) , (6.7) 



and 



Thus, 



5o = -[l + T(C2fe-i,a;,2/)] 



b.^x^'^-^y 



1 r 



Afe.oj — 2 



60 ± \ hl - 46i 



(6.9) 



(6.10) 



3.11) 



By a generalization of Eq. (2.15) in jla | from chromatic 
polynomials to the full Tutte polynomial, it follows that 



T{Sk,m,x,y) 



(aoAfc,o4 + ai) 
(Afe,o,i ^ Afe^o,2) 



(Afe,o,i)'^ 



(aoA/c,o,: 



ai) 



(Afc,o,: 



^fc.O.lJ 



(Afe,0,2)'^ 



(6.12) 



Note that T{Sk,m,x,y) is symmetric under the inter- 
change Afe^o.i ^ Afc_o,2- It is straightforward, using Eq. 
(12. 4p . to re-express these results in terms of the Potts 
model partition function Z(5fc^m,q, u); for brevity, we 
omit the explicit results. 



VII. LONGITUDINAL HOMEOMORPHIC 
INFLATIONS OF CYCLIC AND MOBIUS 
LADDER GRAPHS 

In this section we present an exact calculation of the 
Tutte polynomial for longitudinal homeomorphic infla- 
tions of cyclic and Mobius ladder graphs. We start with 
a cyclic ladder strip of length m squares and add fc — 2 
degree-2 vertices, with fc > 3, to each longitudinal edge. 
This yields the strip graph with longitudinal homeomor- 
phic inflation that we denote Lk^m- The corresponding 
Mobius strip Mk,m is obtained by cutting the cyclic graph 
at any transverse edge and reattaching the ends after a 
vertical twist. The graphs Lk^m and Mk.m each have the 
number of vertices 



n{Lk,m) = n{Mk^,n) = 2(fc - 
and the number of edges 



l)m 



e{Lk,m) = e{Mk^rn) = (2fc - l)m 



(7.1) 



(7.2) 



For a given m and for k = 2 these are the original cyclic 
and Mobius ladder strip graphs; the case fc = 3 is the 
first homeomorphic infiation of these respective graphs, 
and so forth for higher values of fc. 



For the cyclic strip graph Lk^rn we calculate 

2 nr(2,d) 

r(Lfc,„,x,y) = — ^^c^"^) J2 (AfcA.r (7.3) 



d=0 



where c^") = 1, c^^) = q-1, c^^) = q^-3q+l^ nT(2,0) = 
2, 71^(2, 1) = 3, and nr(2 ,2) = 1. The Xk,o,j for j = 1, 2 
were given above in Eq. (|6.1ip . For the others, we find 



fc-i 



where 



Uk 



Afcii = X 



-(life ± ) , j = 2, 3 , 



(7.4) 
(7.5) 



and 



and finally. 



k-2 

3=0 

^k ^k-i _ 2 

X — I 



ru^ul- Ax^ , 



Afc,2,l — 1 



(7.6) 



(7.7) 



Although our result (|7.3p (and ()7.9p below) are formally 
rational functions in a;, one easily verifies that the pref- 
actor 1/(2:: — 1) divides the expression to its right, so that 



9 



T{Lk,rmX,y) and T{Mk^m,x,y) are polynomials in x as 
well as y, as guaranteed by Eq. (|2.2I) . It is again straight- 
forward, using Eq. (j2.4[) . to re-express these results in 
terms of the Potts model partition function Z{Sk,m, Q, v)- 
For the Mobius strip Mk,m, we find that the A's are 
the same as for the cyclic strip Lk,m, and the pattern 
of changes in the coefficients is the same as was shown 
in [lb, 30, 31] for lattice strips without homeomorphic 
expansion, so that 



T{Mk,m,x,y) = 



X - 1 



J=2 



(7.9) 



From these general expressions, one can specialize to the 
chromatic polynomials and the flow polynomials. Via Eq. 
(I4.12p . one readily checks that the results for the chro- 
matic polynomials agree with those that we calculated 
before with S.-H. Tsai in Ref. The flow polynomials 
are unaffected by homeomorphic expansion, and hence 
coincide with those that we calculated with S.-C. Chang 
in Ref. 1331. 



VIII. VALUATIONS OF TUTTE POLYNOMIALS 
FOR HOMEOMORPHIC EXPANSIONS OF 
CYCLIC AND MOBIUS LADDER GRAPHS 

Special valuations of the Tutte polynomial yield several 
quantities of graph-theoretic interest. In this section we 
calculate these for the longitudinal homeomorphic expan- 
sions of cyclic and Mobius ladder graphs. We first recall 
some definitions. A tree graph is a connected graph with 
no circuits (cycles). A spanning tree of a graph G is a 
spanning subgraph of G that is also a tree. A spanning 
forest of a graph G is a spanning subgraph of G that may 
consist of more than one connected component but con- 
tains no circuits. The special valuations of interest here 
are (i) r(G, 1,1) = Nst{G), the number of spanning 
trees (ST) of G; (ii) T{G, 2, 1) = Nsf{G) the number of 
spanning forests {SF) of G; (iii) T{G, 1,2) = Ncssg{G), 
the number of connected spanning subgraphs (GSSG) of 
G; and (iv) r(G,2,2) = Nssg{G) = 2<^\ the number 
of spanning subgraphs (SSG) of G. The last of these 
quantities is determined directly from Eq. (|7.2p as 



SSG 



(8.1) 



where we introduce the notation Gk,m to stand for either 

Lk,Tn or Mk,m- 

We calculate 



A 



= 2 



= 2 



2/0-1 _|_ / 22('^-i) 



1 



1/2 



where the ± sign applies for j = 1,2, respectively. 



= 2 



k-l 



•2) 



(8.3) 



and 



Afc 1, 



^=2 



3-2''-^ -1± 



(3-2''-^ -1)2-2'=+^ 



1/2- 



(8.4) 



where the ± sign applies for j = 2,3, respectively. Next, 
define the notation 



VG = VGk. 



-1 if Gk.m — Mk,m 



(8.5) 



In terms of these quantities, we have, for Gk.m — Lk.-, 
or Mk,m, 

NsF{Gk,m) — (Afe^o,!)™ + {^k,0,2)'" 



-Vg 



l-(Afe,i,i)'" - (Afe,i.2)" + (Afe,i.3)^ 



(8.6) 



where in Eq. (|8.6p the Afc.^j 's are evaluated at x = 2, y = 
1 as in Eqs. (I8.2p - (|8.4|) . For k = 2, i.e. the original cychc 
and Mobius strips, one checks that eq. (|8.6p reduces to 
eq. (D.14) of our previous work [10]. As an example, for 
the first homeomorphic expansion, k ~ 3, eq. (|8.6p yields 

NsF{G3^m) = 

[4(4 + VTb )]™ + [4(4 - VTd )]'" + t]g{1 - 2'") 

'ii + ^/TosA" . /ii-%/T05 



(8.7) 



For valuations of T{Gk,m, x, y) with x — 1 and Gk,m — 
Lk^m or Mk,m, a useful equality is 



Afe,oj — Afe.ij+i 

x—l 



2k-l + y± 



(2k - \ + yf - Ay 



1/2 



x=i 2 

where the ± sign applies for j — 1, 2, respectively. For the number of spanning trees on these families on Gk,m = -^fc,i 
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or M, 



k.m: 



we find 



NsT{Gk,m) = m{k - 1) 



Vg 



(8.9) 



For fc = 2, this reduces to Eq. (D.13) of 
Next, we introduce tlie shortliand notation 



A± = i(2fc + l±\/4P+4A;-7 



and 



a± 



k-1 



k± 



(2fc2 + fc - 4) 
V4fc2 + 4fc - 7 



(8.10) 



(8.11) 



In terms of these quantities, we find, for the number of 
connected spanning subgraphs. 



N, 



CSSG 



(Ah 



l-k+{X+y"-'a+ + iX-) 



\ rn — 1 , 



.12) 



and 



NcssaiMk.rn) = NcssGiLk,^) + l + ^k-l)m . (8.13) 

For fc = 2, this reduces to Eq. (D.15) of pf. 

These graphical quantities grow exponentially as a 
function of the strip length m and hence also n. For 
each quantity Ng (i.e., Nst, Nsf, Ncssg, and Nssg), 
one can thus define a growth constant 



(8.14) 



In Ref. [lO'l we used an equivalent quantity to describe 
the exponential growth, viz.. 



.15) 



We find that for a given set of subgraphs G{G), these 
growth constants are the same for the m — > cxd limits of 
the Sk.mi Lk^m, and Mk^m families of strip graphs, i.e.. 



'^Ng{{Sk}) - ^Ws({ifc}) = '^Ng({Mk)} 



.16) 



so we will just label them by {Lk}- For each type of 
subgraph Q, the growth constant is determined by the 
dominant \k.d,ji which is Afc_o,i, evaluated at the respec- 
tive values of {x,y). We find 



2fc-l 

7/1 — 92(fc-l) 



l+fl-2-2(fe-i) 



1/2 



.17) 



(8.18) 



w 



2fc + l + (4fc2+4fc- 7)1/2 



(8.19) 



TABLE I: Values of 



, , and ui „„ for the m — > oo 



limit of the Gj. , 

Lk,m, or Mk.m- 



SSG' "'S-F' CSSG^ ^ST 

family of graphs, where G^.m denotes Sk ; 



k 


^SSG 








2 


2.8284271 


2.7320508 


2.1357792 


1.9318517 


3 


2.3784142 


2.3689169 


1.6089554 


1.5537740 


4 


2.2449241 


2.2434544 


1.4360948 


1.4104463 


5 


2.1810155 


2.1807485 


1.3466467 


1.3318300 


6 


2.1435469 


2.1434946 


1.2908358 


1.2811894 


7 


2.1189262 


2.1189154 


1.2522226 


1.2454451 


8 


2.1015133 


2.1015110 


1.2236924 


1.2186716 


9 


2.0885476 


2.0885471 


1.2016292 


1.1977615 


10 


2.0785185 


2.0785183 


1.1839857 


1.1809157 


oo 


2 


2 


1 


1 



and 



= fc 



(fc + v/fc2 - i)- 



(8.20) 



For fc = 2, i.e., the original ladder strip without home- 
omorphic expansion, Eqs. (|8T7l) . (jSTSl) . dHUH), and 
(11201), together with reduce to Eqs. (D.24), 

(D.22), (D.23), and (D.21) of our previous paper, Ref. 
[10|. Some numerical values of these growth constants 
obtained from the analytic results given above are dis- 
played in Table HI 



IX. TUTTE POLYNOMIALS OF HAMMOCK 
GRAPHS 

A. General Calculation 

A hammock graph Hk^r is defined as follows: start 
with two vertices connected by r edges (where r denotes 
"rope" or, more abstractly, "route"). Now add fc — 2 
degree- 2 vertices to each of these edges, with fc > 3, so 
that on any rope there are a total of fc vertices, including 
the two end- vertices. Thus, Hk^r with fc > 3 is a uniform 
£ = k — 2 fold vertex (homeomorphic) inflation of the 
original graph, -ff2,r with r edges connecting the two end 
vertices. The number of vertices, edges, and (linearly 
independent) cycles in this graph are 



n{Hk,r) = 2 + (fc - 2)r , 
e{Hk,r) = (fc - l)r , 



(9.1) 
(9.2) 
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and 



(9.3) 



(satisfying the general relation c(G) — e(G) + kc{G) — 
n{G)). Hence, 



Keff{Hk,r) 



2r{k- 1) 
2 + {k~- 2)r 



(9.4) 



An interesting feature of the Hk^r graphs is that if r — ?> 
oo for fixed fc, the combination of tne two end vertices, 
each with k = r and the (fc — 2)r interior vertices on the 
"ropes", each with k = 2, yields the result 



lim Keff{Hk,r) = 2 



1 



(9.5) 



Although Hk^r is much simpler than the complex net- 
works encountered in biological and social contexts, the 
property that for large r it has vertices of quite differ- 
ent degrees is also observed in complex networks 's^. 
The girth g{G) of a graph G is the number of edges in a 
minimum-distance circuit in G. For r > 2, the girth of 
Hk,r is 



g(Hk,r) = 2k~2 



(9.6) 



Before presenting our general results for T{Hk^r, x, y), 
we note two special cases. First, for r — 2, the graph 
Hk,2 is just the circuit graph with 2k — 2 vertices: 



k,2 



C2k- 



Hence, 



T{Hk^2,x,y) = T{C2k-2,x,y) , 



(9.7) 



(9.8) 



where T{Cn,x,y) was given in Eq. (|6.6|) . The second 
special case is for k = 2. If G is a planar graph, we 
denote its planar dual as G*. From Eq. (|2.2|) . it follows 
that 



T{G,x,y)^T{G*,y, 



(9.9) 



We can apply this result here, since Hk^r is a planar 
graph. Now 7J2,r is the graph consisting of two vertices 
connected by r edges. This is the planar dual to the 
circuit graph; 

H2.r = (Cr)* . (9.10) 

From Eqs. (H!]), (|9J0l) . and (HH), it follows that 



r(i?2,.,x,2/)=T(G.,y,x) = ^-t^^-^^ . (9.11) 



2/-1 



Proceeding to the general Hk.r graph, we define 



k-2 



Ai/,1 = J2 ■ 



3=0 



x'^-i - 1 

X- 1 



(9.12) 



and 



Ah^2 = T(Gfc_i,x,?/) 



(9.13) 



Using an iterative application of the deletion-contraction 
relation, we calculate the Tutte polynomial of H^.r to be 



T{Hk,r, X, y) = (A_y,i)'' T{C2k-2,x, y) 



(A//,i) 



r-2 



(Aif,2) 



r-2 



A/fa — \h,2 



(9.14) 



r 



B. Chromatic and Flow Polynomials 

Using Eq. (|4.12p . one readily verifies that for the case 
X = I — q, y — 0, the Tutte polynomial (I9.14p yields the 
chromatic polynomial for Hk,r that we calculated in (Eq. 
(3.7) of) Ref. 

One can also discuss other special cases. For the flow 
polynomial F{Hk.r, q), we set x = 0,y = l — q and, with 
eq. (|4.19|) . we have 



F{H,,r,q)^q-'[{q-ir + {q-l){-m 



(9.15) 



dual, then the flow and chromatic polynomials satisfy the 
relation 



FiG,q)=q-'P{G*,q) 



(9.16) 



We observe that F{Hk^r,q) = q^^P{Gr,q), m accord 
with (|4:20l) and the fact that iJa = (G^)* (cf. Eq. 



C. Reliability Polynomial 



This is independent of /c, in accordance with the general 
result (|4.20l) . If G is a planar graph and G* is its planar 



A communication network, such as the internet, can 
be represented by a graph, with the vertices of the graph 
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representing the nodes of the network and the edges of 
the graph representing the communication hnks between 
these nodes. In reahstic networks, both the nodes and 
the hnks between them are imperfect, and fail to oper- 
ate. One common measure of the rehability of the net- 
work is the probability that there is a working commu- 
nications route between any node and any other node. 
This is the all-terminal reliability function. This is com- 
monly modeled by a simplification in which one assumes 
that each node and link are operating with respective 
probabilities Pnode and punk- As probabilities, Pnode and 
Piink lie hi the interval [0,1]. The dependence of the all- 
terminal rehability function Rtot{G,Pnode,PUnk) on Pnode 
is an overah factor of {pnodeY; i.e., Rtot{G,pnode,Piink) = 
(Pnode)"^ R{G , piink) ■ Thus, the difficult part of the cal- 
culation of Rtot{G ^Pnode^Piink) IS thc determination of 
R{G iPiink)- For notational brevity, we set punk = P- 
The function R{G,p) is given by 

RiG,p) = P'^^^ (1 -p)^(G)-e(G) ^ (g_-^7) 

GCG 

where G is a connected spanning subgraph of G. Clearly, 
R{G,p) is a monotonically increasing function of p e 
[0, 1] with the boundary values R{G, 0) = and 
R{G,1) = 1. R{G,p) can be related to a special case 
of the Tutte polynomial, evaluated with x — 1 (guaran- 
teeing that G is a connected spanning subgraph of G) 
and y = 1/(1 — p)- This relation is 

R{G,p) = p"-i(l - 1 ) (9_i8) 

1-p 

As in our previous calculations for lattice strips we 
can thus obtain reliability polynomials as special cases of 
Tutte polynomials. Using our calculation in Eq. (I9.14p of 
T{Hk,r,x,y) with x = 1 and y = 1/(1 -p) in Eq. ((9l8| . 
we have calculated R{Hk^r,p)- For r = 2, we have 

R{Hk,2,p) - R{C2k-2,p) = p''"'b + 2(fc - m-p)] , 

(9.19) 

in accord with Eq. (19. 7p . For k — 2, i.e., the case of 
no homeomorphic expansion, we find, in accord with Eq. 
dmni), that 

R{H2,r,p) = R{{GrY,p) = l-{l-pr . (9.20) 

In general, we find that for fixed p e (0, 1) and fixed r, 
R{Hk^r,p) is a monotonically decreasing function of k. 
This can be interpreted as a consequence of the fact that 
as k increases, the girth g{Hk,r) increases (cf. Eq. (|9.6p ). 
and hence there is a greater likelihood that one of the 
communication links along the minimum-distance path 
and other paths between two nodes is not operating. 

In contrast, for fixed p € (0, 1) and fixed fc > 3, we find 
a variety of behaviors for R{Hk.r,p) as a function of r. 
To illustrate this, we take fc = 3 and compare a few pairs 
of values (r, r') = {r,2). We calculate 

i?(i73,2,p) =p'(4-3p) , (9.21) 



i?(if3.3,p) = p\7p^ - 18p 12) , (9.22) 

RiH3,4,p) = /(4 - 3p)(5p2 _ I2p + 8) . (9.23) 

and 

R{Hz,^,p) =/(3l/-150p^-H280p2-240p-H80) (9.24) 
Hence, for the pairs (r, r') = (r, 2), 

i?(il3,3,p)-i?(i?3,2,p)=p'(l-p)'(7p-4) , (9.25) 

i?(H3,4,p)-i?(i?3,2,p) =p'(l-p)'(4-3p)(5p2-2p-l) , 

(9.26) 

and 

i?(H3,5,p)-i?(i?3,2,p) = p^(l-p)'(31p5-88/+73p3_6p2-5p 

(9.27) 

As is evident from the difference in Eq. (|9.25p . 

RiH3,3,p)> R{H3a,p) if 1>P>^ -0.571, 

(9.28) 

while 

R{H^,z,p) < R{Hza,P) if < p < ^ . (9.29) 
Similarly, Eq. (|9.26p shows that 

R{H3,4,P) > R{H3,2,P) if l>p>i^^~ 0.690, 

(9.30) 

while 

RiH3A,p) < R{H3,2,P) if < p < . (9.31) 

5 

Similarly, i?(i?3,5,p) > i?(i?3,2,p) for 1 > p > 0.8418, 
and i?(iJ3,5,p) < R{H3,2,p) for 0.8418 > p > (where 
the crossover value of p is a root of the quintic in Eq. 
(I9.27P quoted to four significant figures). Note how, for a 
fixed value of fc, the value of p beyond which R^H^^riP) 
is greater than R(H3^2,p) increases with r, from 0.571 
for r = 3 to 0.690 for V = 4 to 0.842 for r = 5 (to three 
figures accuracy), and so forth for higher values of r. We 
observe similar behavior for R{Hk^r,p) — R{Hk,2iP) as a 
function of r for higher values of fc. 
However, we also have 

i?(i?3,4,p)-i?(i?3,3,p) --/(l-p)'(15p'-26p+12) , 

(9.32) 

and 

R{H3.5,p)~RiH3^i,p) = V(l-p)'(-31p3+88p2_88p+32) . 

(9.33) 

These comparisons provide a contrasting type of behav- 
ior, since 

R{Hs,5,p) < RiH^A.p) < i?(i?3.3,p) V p e (0, 1) . 

(9.34) 

and so forth with R{H3,r,p) for larger values of r. Thus, 
for these cases, for all p e (0,1), R{H3^r,p) is a mono- 
tonically decreasing function of r as r increases above 
3. 
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D. Some Graphical Quantities 



and 



We next discuss special valuations of T{Hk.r, x, y) that 
yield quantities of graph-theoretic interest. For x — I, 
the following result is convenient: 



T{Hk^r.. 1, y)^{k- lY-^(2k + y - 3) 
(fc + y-2)''-2_(fc-i)'-2- 



y-1 



ik + y-2y 



In addition to 



we compute 



NsTiHk,r) ^ rik - ly 



(9.35) 
(9.36) 
(9.37) 



NsFiHk,r) = (2'=-i-l)'-2r22('=-i)-l+(r-2)(2'=-i-l) 

(9.38) 



and 



N, 



CSSG 



{Hk,r) 



(k-iy 



(9.39) 



Since these are two-parameter families of graphs, one 
can consider the limits (i) r — >■ (X) with fixed finite fc, 
and (ii) A: — >■ oo with fixed finite r. We calculate the 
growth constants for each of these. The cases of interest 
here are those with k > 3, i.e., those with homeomorphic 
expansion. For these we find that for the limit (i) , 



2 — 



(9.43) 



(9.44) 



(9.45) 



(9.40) 



For the limit (ii), we calculate 

(a) : w^sG = u'sF = 2 

and 

(ii) : = w^ssG = 1 

X. CONCLUSIONS 



In conclusion, in this paper have have derived exact 
relations between the Potts model partition function, or 
equivalently, the Tutte polynomial, for a graph G and 
for a graph G obtained from G by edge or vertex infla- 
tion. An analysis was given of some physical effects of 
uniform ^-fold edge and vertex inflations. We have pre- 
sented exact calculations of the Tutte polynomials for 
free, cyclic, and Mobius ladder families of graphs Sk,m, 
Lk^rm and Mk^m of length to, with i = k — 2 vertex 
inflation on all longitudinal edges. We have given simi- 
lar results for the Tutte polynomial of the family Hk.r of 
hammock graphs. Our present results generalize our pre- 
vious calculations in Refs. 0, flllSl- As one application, 
we have calculated reliability polynomials for the ham- 
mock graphs and analyzed their properties as a function 
of k and r. In addition, we have used our calculations to 
compute the number of spanning trees, spanning forests, 
and connected spanning subgraphs on these families and 
to determine their asymptotic behavior as the number of 
vertices goes to infinity. 



(z) : w,^ = ik-l)^ , (9.41) 
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